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Sub-pixel resolution with color X-ray camera SLcam® 

Stanistaw H. Nowak,* a Aniouar Bjeoumikhov/^ Johannes von Borany, c Josef Buchriegler/ Frans 
Munnik, c Marko Petrie/ Martin Radtke/ Axel D. Renno/ Uwe Reinholz/ Oliver Scharf, a and 
Reiner Wed ell" 

The color X-ray camera SLcam® is a full-field, single photon detector providing scanning free, energy and spatially resolved 
X-ray imaging. Spatial resolution is achieved with the use of polycapillary optics guiding X-ray photons from small regions on 
a sample to distinct energy dispersive pixels on a charged-coupled device detector. Applying sub-pixel resolution, signals from 
individual capillary channels can be distinguished. Accordingly the SLcam® spatial resolution can be released from pixel size 
being confined rather to a diameter of individual poly capillary channels. In this work a new approach to sub-pixel resolution 
algorithm comprising photon events also from the pixel centers is proposed. The details of the employed numerical method and 
several sub-pixel resolution examples are presented and discussed. 


1 Introduction 

SLcam® is a high quantum efficiency and throughput color X- 
ray camera system 1 3 designed for divergent X-ray radiation. 
It allows detection of single photons with both energy and spa¬ 
tial resolution. SLcam® combines a pn-junction Charged- 
Coupled Device (pnCCD) 4 with polycapillary optics] 5 The 
latter can be regarded as a bunch of independent X-ray chan¬ 
nels guiding X-ray photons from small regions on a sample to 
corresponding pixels on pnCCD similarly to the way as fiber 
optics guide light. 

SLcam® employs pnCCD with 48 fim pixel size. Keeping 
the pixel size fixed spatial resolution of SLcam® can be se¬ 
riously improved by the use of conically shaped magnifying 
optics.B 3 This type of optics is available with magnification 
factors up to 10:1, allowing representation of 4.8 x4.8 gm 2 
area in a single pixel. Theoretically, according to Nyquist- 
Shannon sampling theorem,® such a system can correctly re¬ 
solve details down to ~ 9 gm. 

Currently used SLcam® optics are optimized for the 
pnCCD pixel dimension. The capillary exit diameter is 
adapted in such a way that the spot size from an individual 
channel on the detector is approximately equal to the pixel 
size. However, current technology allows fabrication of poly- 
capillaries with single channel diameter in the range of a mi¬ 
cron or even below giving room for further improvement of 
resolution. Currently pixel size is the limiting factor for spa¬ 
tial resolution of SLcam®. 


a IfG - Institute for Scientific Instruments GmbH, Berlin, Germany. E-mail: 
nowak@ ifg-adlershof.de 

b Institute for Computer Science and Problems of Regional Management, 
Kabardino-Balkaria, Russia. 

c Helmholtz-Zentrum Dresden-Rossendorf Dresden, Germany. 
d J. Stefan Institute, Ljubljana, Slovenia. 

e BAM Federal Institute for Material Research and Testing, Berlin, Germany, 
f Helmholtz-Zentrum Dresden-Rossendorf, Helmholtz Institute Freiberg for 
Resource Technology, Freiberg, Germany. 

8 IAP Institute for Applied Photonics e.V., Berlin, Germany. 


With the use of a sub-pixel algorithm the dominant role of 
pixel size can be released. This algorithm divides the signal 
assigned to each physical pixel over a number of virtual sub¬ 
pixels. With such an approach further downscaling of poly¬ 
capillary channels can be practically used for SLcam® lateral 
resolution improvement. 

Images with sub-pixel resolution are achievable due to spe¬ 
cific physics of X-ray photon with Charged-Coupled Device 
(CCD) interaction that guides to creation of a so-called elec¬ 
tron cloud. A non zero area of the cloud leads to charge depo¬ 
sition in pixels nearest to the photon hit. With a correct recon¬ 
struction of the footprint of a single photon event the photon 
hit position can be estimated with a much higher precision 
than the pixel size. 

The sub-pixel resolution algorithm was first applied to Ad¬ 
vanced CCD Imaging Spectrometer (ACIS) installed at Chan¬ 
dra X-Ray Observatory (CXO). 7 12 Some years later similar 
technique was adapted to pnCCD. 13 ^Unfortunately, the pro¬ 
posed routines are neglecting a majority of events analyzing 
only so-called ’’comer events”, i.e., the photon hits that reach 
the CCD in proximity to a pixel corner. As a result the sub¬ 
pixel image exhibit a grid pattern with the intensity drops in 
pixel centers. El We propose a modified approach taking all the 
photon hits and distributing them properly over the sub-pixel 
pattern. 

2 Method 

2.1 Electron cloud 

An X-ray photon absorbed in silicon generates a number of 
electron-hole pairs amounting to E /Wsi - the photon energy E 
divided by the formation energy of a single electron-hole pair 
W Si = 3.6 eV Pin a fully depleted layer the carriers are sepa¬ 
rated in the vertical electric field and diffuse laterally produc¬ 
ing two charge clouds with opposite signs. The holes are col- 







lected at the large area cathode. The electrons are transfered 
to the pixelated anode where they are split over individual pix¬ 
els. During the readout these charge packets are sequentially 
transferred to charge amplifiers and counted. 

The carrier density of a single photon electron cloud is fre¬ 
quentl y app roximated by a two dimensional Gaussian distri¬ 
bution 11 3 18 : 
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where g is the total charge produced, o x and o y are the v and 
y widths of the cloud, and (xoDo) 1S P°i nt of a photon hit. 
For a pixel with coordinates (z, j) the resulting charge intensity 
7j can be calculated as an integral over a square corresponding 
to a pixel area: 
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Note that 7j is a product of two separable functions depending 
only on one variable x or y. 

Though this is a simplified model neglecting the energy de¬ 
pendence, charge quantization, etc ., it shows the fundamental 
fact making the sub-pixel resolution algorithm possible: on a 
pixelated plane a single photon electron cloud is normally split 
over several pixels. For every photon hit, in addition to pixel 
position, extra information in a form of pixel intensity distri¬ 
bution is gathered. This additional information can be used to 
determined the center of impact of the photon with sub-pixel 
accuracy. 


2.2 Intensity ratios 


A good measure of a single photon electron cloud distribution 
are the intensity ratios defined as follows: 
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where I\ and I 2 correspond to charge gathered in two different 
areas of a CCD. Possible values of d(I\ .I 2 ) ranges from —1 to 
1; 0 is reached when I\ = I 2 . 

It can be shown that for any electron cloud distribution for 
which x and y components can be separated ( e.g ., Gaussian 
distribution) the pixel intensity ratio in x direction is indepen¬ 
dent from y coordinate and vice versa : 
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Note also that, provided that for each photon hit the shape of 
the electron cloud is uniform, the intensity ratio should change 
monotonically with v and y coordinates of the hit position. 

A very elegant way to estimate the single photon electron 
cloud distribution was presented in Ref [19] The method relies 
on analysis of a histogram of measured intensity ratios and 
assumes only a uniform distribution of photon hits over the 
CCD plane and a Gaussian shape of the charge cloud. The 
method was used to reveal the size of the cloud. We will use 
similar methodology to calculate sub-pixel coordinates. 

In order to calculate the intensity ratio histogram 2x2 pixel 
boxes centered at a pixel corner (x^yj) nearest to the point of 
a photon hit (voDo) are analyzed; intensity ratios in, respec¬ 
tively, v and y directions are calculated as follows: 

dx = d{Pr\+pr\ 1)^1)) (6) 

d y = diVrX+I)^, iy'+I))- (7) 

The intensity ratios in x or y direction computed for all the 
single photon hits are combined in one histogram. 

The histogram normalized to occupy a unitary area can be 
regarded as an estimation of the probability density function. 
Accordingly, a normalized intensity ratios histogram rates the 
probability of a single photon hit to create a given intensity 
ratio value. 

In Fig.[l]an example of normalized intensity ratio histogram 
is plotted. The figure shows data computed for Au L pho¬ 
tons from the measurement of Au bar pattern presented fur¬ 
ther in Section [3] As can be seen most of the photon hits 
create electron clouds with the intensity ratio in proximity to 



Fig. 1 An example intensity ratio histogram of Au L photons 
generated from a real measurement (Au bar pattern) with 
formula The histogram is normalized to occupy a unitary area 
and can be identified with the probability density function. In the 
inset the histogram is compared to Gaussian simulation (red line). 




















— 1 or 1. Intensity ratios close to 0 are the least probable. 
Discussed Gaussian model fits well to most values of the his¬ 
togram; though for | d x \ ~ 1 (and \d y \ « 1) the simulation fails 
showing deficiency of the model. 


2.3 Sub-pixel coordinates 

Our method consist in converting the d x and d y intensity ra¬ 
tios of a given photon hit to vo and yo coordinates. In our 
approach we do not model the shape of the single photon elec¬ 
tron cloud; we assume only separation of v and y components 
of the charge distribution. We also assume a uniform distri¬ 
bution of photon hits over the CCD plane. Accordingly the 
relation between the intensity ratio and single photon hit posi¬ 
tion is not simulated, but measured. 

Lets set (v^yy) - the pixel corner nearest to the point of a 
photon hit - to the origin, and the number of photon events 
to N. The possible values for vo inside a 2x2 pixel box span 
within (— f j §)» w here P is the pixel dimension. Due to uni¬ 
form distribution of photon hit positions over the CCD plane 
the probability of finding a photon hit with x coordinate below 
vo has a linear form: 



Fig. 2 Relative photon hit position vo calculated as a function of 
intensity ratio d x according to equation fj~3j ). Likewise in Fig.[l]the 
data for Au L photons from Au bar pattern are illustrated. 
Secondary vertical axis shows cumulative probability distribution 
Pd x . In the inset a 2 x 2 pixel box with 3x3 sub-pixel division is 
shown. Different levels of gray mark - in order from darkest to 
brightest - central, left and right sub-pixels. 
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here n() stands for the number of elements of a set. We can 
also compute cumulative probability functions of d x and d y : 
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P dx {d x ) and P dx (d x ) are equivalent to normalized cumulative 
histograms of the measured intensity ratios. 


If we note that d x is monotonically increasing with vo then 
we get: 

n(d x : d x < d x (x 0 )) =n(x: v<v 0 ). (11) 

This directly shows that P dx can be identified with P x \ 


P dx (d x (x 0 )) = P x (x 0 ). (12) 

Finally an explicit relation for vo can be given: 

xo = p (pdMx) - 0 ■ (1 3 ) 

An analogical relation is also valid for yo and d y . 

The relationship between vq, d x and P dx is shown in Fig. [2] 


2.4 Division to sub-pixels 


A strong point of our approach is that it requires only a good 
statistics of single photon hits ensuring high enough number 
of events N to minimize the statistic error of generated P dx and 
P dy . In particular the method is independent from the shape 
and creation process of the electron cloud and can deal with 
any, also asymmetric, 9 charge distributions. The cloud cre¬ 
ation process and the resulting charge distribution is a charac¬ 
teristics of a given CCD and should not change in time. There¬ 
fore each measurement could be analyzed with the use of the 
same, previously obtained pixel intensity ratio probabilities. 

In an ideal case of uniformly shaped electron clouds and 
a hypothetical noise-free CCD, equation f]~3] ) should give a 
strict position of a single photon hit. In reality the electron 
cloud is shaped randomly and the intensity signal from each 
pixel is given with an error. These two factors contradict the 
assumptions on monotonic relation between d x and vo and on 
separation of v and y components. As a result the accuracy of 
sub-pixel coordinates is limited. 

We can estimate an error on d(I\ ,h) with the following re¬ 
lation: 
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where 81 is an error on a charge gathered in I\ or 1 2 . Note 
that 8d(l\ .I 2 ) has its maximum for I\ = 0 or I 2 = 0 which 
corresponds to | d(I\ ,h) | = 1 or |vo | = §. This is exactly where 
P dx exhibits the largest increase and, respectively, where the 






























error on P^ x is the biggest. 

As a result the photon hit position is most inaccurate in 
proximity to the pixel center |jco| = f . In order to minimize 
that effect the sub-pixel division should comprise that higher 
uncertainty central region in a single sub-pixel. An exemplary 
division to 3 sub-pixels in v direction is shown in Fig. [2] 

In conclusion, the assignment of sub-pixel coordinates to a 
given photon event consist only of two steps: (i) calculation 
of d x and d y with formulas ([6| and (0, and (ii) comparison 
of obtained values with sub-pixel borders expressed in pixel 
ratio coordinates. As a result the algorithm is fast and can be 
applied on-line. 

3 Experimental 

To examine the performance of the sub-pixel algorithm two 
test samples were measured. The first one contains several 
microns of Au on a Si sup por t arranged in bar-like patterns and 
a uniform reference layer.^ 0 The bar pattern consists of 10 Au 
lines with the width and the spacing between them decreasing 
from 10 /im to 1 /u m in a 1 /im steps. The resolution limit can 
be easily found as the width of the narrowest recognizable line. 
The second structure was produced by the Fraunhofer-Institut 
fur Zuverlassigkeit und Mikrointegration (IZM) in Berlin and 
consists of 3 /im thick and 30 fim wide Cu stripes deposited 
on Si wafer with 200 nm TiW adhesion layer. The stripes 
are aligned in parallel and distributed in groups with constant 
spacing of 30, 50 and 90 /im. Sketches of both structures are 
illustrated in Fig. [3] 

Imaging was performed with two types of polycapillary 
lenses - a high-resolution conical optics with a 8:1 magnifi- 
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Fig. 3 Sketch of Au bar pattern (left) and Cu stripes (right). 
Structures were measured with, respectively, focusing 8:1 and 
parallel 1:1 optics; sizes of corresponding pixels (side length of, 
respectively, 6 /rm and 48 jUm) and capillary channels (entrance 
diameter of, respectively, 2 fim and 24 /im) are shown for 
comparison. 


cation that was employed for the Au pattern, and a parallel 1:1 
optics used for the Cu stripes structure. The magnification op¬ 
tics has a single channel entrance diameter of di n = 2 /im. The 
exit diameter of a channel span to d out = 16 /im which is one 
third of the pixel size (p = 48 /im). A single channel diameter 
of the 1:1 optics amounts to d = 24 /im which is half of the 
pixel size. 

In order to ensure high enough photon count rate, mea¬ 
surements were performed at synchrotron radiation facility 
with the beam provided by the BAM line at BESSY 11“ 
and at newly developed PIXE beam line (HS-PIXE) at 
Ion Beam Center at Helmholtz-Zentrum Dresden-Rossendorf 
(HZDR). 22 Detailed descriptions of two setups can be found 
in Ref. GO and [22] 

For imaging purpose Au L and Cu K lines intensity dis¬ 
tributions were assessed. Valid photon events were selected 
with an algorithm described elsewhere. 3 The algorithm rejects 
the noise such as cosmic rays and accepts only photon events 
having an appropriate arrangement of pixels above the noise 
threshold. 

The position of the photon hit is first estimated by the 
weighted position of the pixel with the highest intensity and 
its nearest neighbors. Subsequently a 2x2 pixel box around 
a pixel comer closest to the estimated position is selected and 
the sub-pixel algorithm is applied. The probability functions 
of pixel intensity ratios in v and y directions were obtained di¬ 
rectly from the measurements. Probability density and cumu¬ 
lative probability functions of the intensity ratios in v direction 
for Au L photons are presented in Fig. [I] an d0 

4 Results 

4.1 Au bar patter 

In Fig. [4] the image of Au bar pattern is presented in normal 
resolution and with pixels divided to 2x2 and 3x3 sub-pixels. 
The normal resolution image exhibits strong pixelation intro¬ 
ducing artificial jerking of the Au bars and decreasing the con¬ 
trast between structure and background. Only the 4 thickest 
stripes can be distinguished in the original image; the rest is 
blurred. With the division into 2x2 sub-pixel another Au line 
is resolved. The 3x3 subpixel division, however, does not 
introduce any visible improvement. In this regime a single 
sub-pixel correspond to 2 /im distance on the sample which is 
exactly the entrance channel diameter of the polycapillary op¬ 
tics (2 /im). In addition the small effective area of a sub-pixel 
significantly lowers the count rate and increases the noise con¬ 
tribution. 

In Fig. [5] three intensity profiles along a single (sub-)pixel 
row spanning along the whole Au bar structure are depicted. 
A profile intensity drop with pixel division is clearly visi¬ 
ble. This is an effect of the sub-pixel area decrease. For a 




























Fig. 4 Normal resolution (left), 2x2 (center) and 3x3 (right) sub-pixel division images of the Au bar pattern. The image was acquired at the 
BAM line at BESSY II with an 8:1 magnification lens. The grid shows the positions of the real pixels, each corresponding to 6x6 gm 2 area on 
the sample surface. Red flashes indicate the (sub-)pixel row chosen for comparison in Fig. [5] 
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Fig. 5 Pixel intensity profiles along one sub-pixel row indicated in 
Fig .[4] with red rims. The profiles are plotted for normal resolution 
(black line), 2x2 (blue line), and 3x3 (red line) sub-pixel division. 
Estimated positions of the Au bar centers are indicated with black 
vertical lines. On top the actual Au stripe widths are given. 
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Fig. 6 Intensity profile across Au bar structures calculated for 
normal resolution (black line) and 3x3 sub-pixel division (red line) 
images. Estimated positions of the Au bar centers are indicated with 
black vertical lines. The intensity profile was constructed from data 
sampled at regular intervals across the Au bars. To increase the 
signal to noise ratio the intensity was integrated over the whole 
length of Au bars. 


























































































































































pixel divided to nxn sub-pixels the intensity is decreased n 2 
times. However, the increased sampling frequency evidently 
increases the contrast. For the profile with no sub-pixel di¬ 
vision only the 4 thickest Au bars can be definitely differen¬ 
tiated. The intensity valley between 6 gm and 5 /im bars is 
not visible. In case of 2x2 pixel division the 6 /im bar can be 
clearly distinguished. For 3x3 sub-pixel division, due to the 
considerable noise content, there is no significant amendment. 

It should be noted that due to the Nyquist-Shannon sam¬ 
pling theorem,® details smaller than twice the sampling dis¬ 
tance cannot be correctly distinguished. The theorem says 
that making n samples over a certain distance only a signal 
comprising less than n /2 elements can be rendered properly. 
Higher number of elements will lead to signal aliasing result¬ 
ing in distortions and artifacts. Of course, with the sub-pixel 
division, this constrain is relaxed. 

As a single pixel in a 8:1 magnified image represents a 
square with a side length of 6 /im the camera in standard reso¬ 
lution cannot correctly represent features smaller than 12 /im. 
Still, as can be seen in the presented example, even without 
sub-pixel resolution much smaller lines can be distinguished; 
though with some alterations. For instance in Fig. [5] the nor¬ 
mal resolution profile shows the 8 /im bar to be thicker than 
the 9 /im one. 

In Fig. [6] two additional intensity profiles are shown. The 
profiles represent the cumulative intensity across the lines par¬ 
allel to Au bars sampled at regular intervals. The sampling 
interval was identical for both curves corresponding to normal 
resolution and 3x3 sub-pixel division; therefore the relative 
intensity of two profiles does not differ. For the sake of leg¬ 
ibility the profile corresponding to 2x2 sub-pixel division is 
not present. 

The tilt of the Au bar structure with respect to the pixel lines 
eliminates the effects due to pixelation. As a consequence the 
6 /im Au bar can be resolved even in case of normal reso¬ 
lution. Also here the contrast improves when sub-pixel res¬ 
olution is applied. In addition the enlarged integration area 
increases the statistics leveling the noise contribution. As a re¬ 
sult the sub-pixel resolution profile allows distinction of 5 /im 
or even 4 /im bars. 

4.2 Cu stripes 

The consequences of too sparse sampling are clearly visible in 
the image of Cu stripes (see Fig. [7]). The image represents a 
structure with three groups of Cu stripes with dimensions be¬ 
low the Nyquist-Shannon resolution limit. Here the structure 
is aligned with the pixels and the effect of aliasing is clearly 
visible. Normal resolution and 3x3 sub-pixel division images 
are presented for comparison. 

The left most group is too confined and individual stripes 
cannot be distinguished even when sub-pixel resolution is ap- 



Fig. 7 Normal resolution (left) and 3x3 sub-pixel resolution (right) 
images of the Cu stripes pattern. The image was acquired at 
HS-PIXE beam line at HZDR with a parallel 1:1 optics. The grid 
shows the positions of the real pixels, each corresponding to 
48 x48 /im 2 area on the sample surface. On top the cumulative 
intensity along each (sub-)pixel column is presented. The dotted 
guide lines indicate the contrast levels for the sub-pixel resolution 
image. 

plied. The structure has stripes with dimensions and spacing 
of 30 /im which is too close to the single capillary channel 
diameter (24 /im). 

The stripes in the central group are distributed with inter¬ 
val of 90 /im, which is comparable to the Nyquist-Shannon 
resolution limit, i.e., double the pixel size (p = 48 /im). In 
the standard image the stripes are merely distinguishable and 
confined to a single pixel line. The sub-pixel resolution image 
shows a much better representation of the structure. 

Finally, the group on the right has an interval large enough 
to fulfill the Nyquist-Shannon theorem, but individual stripes 
are still below the limit. Therefore the structure is clearly re¬ 
solved, however the stripe’s widths and positions are not ren¬ 
dered well giving a false impression of a patchy structure. 
The sub-pixel resolution image correctly presents uniform, 
equidistant stripes. 

For both central and right most group of stripes there is a 
noticeable increase of the contrast level when sub-pixel reso¬ 
lution is used. However, the contrast of the central group is 
much smaller. This is not surprising as the dimensions in the 
central group are much closer to the size of a single polycap¬ 
illary channel which is an ultimate resolution limit. 

4.3 Au layer 

In order to test the limits of sub-pixel resolution algorithm 
an additional assay was performed on a reference Au layer. 
A small portion of the image was compared for the case of 
normal pixel resolution and increasing sub-pixel division (see 
Fig.§. As can be seen very dense divisions lead to creation 














































Fig. 8 Normal resolution and 3 x 3, 5x5, and 9x9 sub-pixel division images of a Au layer. The image was acquired at the BAM line at 
BESSY II with a 8:1 magnification lens. Artificial intensity structure following the pixel center position is strongly indicated for the sub-pixel 
resolution image. The grid shows the real pixel positions. 


of intensity artifacts following the pixel center positions. 

This effect was already discussed in Section 2 A and is a 
result of the uncertainty of a position of a photon hit reach¬ 
ing the middle area of a pixel. Nevertheless, pixel divisions 
up to 5 x 5 sub-pixels are usually not so much affected by this 
inconvenience. It should be noted that due to a nonzero diver¬ 
gence of transmitted photons and a nonzero optics-detector 
distance in SLcam® the footprint of a single polycapillary 
channel on the detector cannot decrease below a dozen of mi¬ 
crons. Thus, keeping in mind that the dimension of a pnCCD 
pixel is 48 jum, a pixel division into 5x5 sub-pixels should be 
sufficient to reach the optimal limits of resolution. 


5 Summary 

We propose a modified approach to sub-pixel resolution algo¬ 
rithm taking into account all the photon events occurring on 
a CCD plane. The sub-pixel position of the photon hit is as¬ 
sessed from the pixel footprint of the generated electron cloud. 
The calculations are performed based on the pixel intensity ra¬ 
tios of a 2x2 pixel box holding photo electrons. The method 
is independent from the actual shape and creation process of 
the charge cloud and does not reject the photon events from 
the pixel center. The algorithm is fast and can be employed 
on-line. 

The sub-pixel resolution was applied to several test struc¬ 
tures. A notable enhancement in quality of the acquired im¬ 
ages comprising contrast and resolution improvement, as well 
as elimination of aliasing due to pixelation was demonstrated. 
For images acquired with 8:1 magnifying optics a resolution 
limit of 5 gm was assessed. Due to an electronic noise and 
random variations of electron cloud shape the sub-pixel co¬ 
ordinates cannot be given with unlimited precision. It was 
shown that in case of 48x48 jum 2 pixel pnCCD the 5x5 sub¬ 


pixel division is a realistic limit of the method. 
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